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Abstract 

We consider a system of 2D interacting fermions with a flat Fermi 
surface. The apparent conflict between Luttinger and non Luttinger 
liquid behavior found through different approximations is resolved by 
showing the existence of a line of non trivial fixed points, for the 
RG flow, corresponding to Luttinger liquid behavior; the presence of 
marginally relevant operators can cause flow away from the fixed point. 
The analysis is non-perturbative and based on the implementation, 
at each RG iteration, of Ward Identities obtained from local phase 
transformations depending on the Fermi surface side, implying the 
partial vanishing of the Beta function. 

Pacs numbers: 71.10.Hf,71.10.Fd 

1 Introduction 

The properties of the 2D interacting fermions are still largely unknown, de- 
spite the tremendous effort devoted to their understanding in the last years. 
One of the most debated questions is on the possible existence of a Luttinger 
liquid phase, first suggested by Anderson [1] as an explanation of some prop- 
erties of high Tc superconductors, as observed also in recent experiments, see 
e.g. [2] . 

It has been proved, in the case of symmetric, smooth and convex Fermi 
surfaces (like in the Jellium model [3] or in the Hubbard model in the non 
half filled case [4]), that the wave function renormalization Z is essentially 
temperature independent up to exponentially small temperatures. As in a 
Luttinger liquid one expects instead a logarithmic behavior in this regime, 
i.e. Z ~ 1 -|- 0([/^log/3), such results rule out for sure the possibility of 
Luttinger liquid behavior. 
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On the contrary the presence in the Fermi surface of flat regions can 
produce non Fermi hquid behavior. The simplest model with a flat Fermi 
surface is the 2D Hubbard model at half filling, in which the Fermi surface 
is a square. It was proved in [5] that the wave function renormalization is 

up to exponentially small temperatures; the presence 
of the log^ /3 is a consequence of the Van Hove singularities, related to the 
fact that the Fermi velocity is vanishing at the corners of the squared Fermi 
surface, and implies that also such a model does not show Luttinger liquid 
behavior. 

It is important to stress that the results in [3], [4], [5] are rigorous as they 
are based on expansions which are convergent provided that the temperature 
is not too low, the finite temperature acting as an infrared cut-offs; how- 
ever such expansions cannot give any information on the zero temperature 
properties. 

A lot of attention has been devoted in recent years to the zero temper- 
ature properties of Fermi surfaces with flat regions and no corners, which 
share some features with the Fermi surfaces of some cuprates as seen in 
photoemission experiments. Parquet methods results [6] and perturbative 
Renormalization Group (RG) analysis [7] truncated at one loop indicate that, 
for repulsive interactions, there is no indication of a Luttinger liquid phase 
at zero temperature; the effective couplings flow toward a strong coupling 
regime related to the onset of d-wave superconductivity. In a more recent 
RG analysis truncated at 2 loops [8] one still gets a flow to strong couphng, 
but in some intermediate region some indication of Luttinger liquid behavior 
is found. 

Apparently conflicting results are found by applying bosonization: in [9] , 
[10] a model of electrons on a square Fermi surface was mapped in a collec- 
tion of fermions on coupled chains, and it is found that the correlations at 
zero temperature in momentum space arc similar to the one of the Luttinger 
model. A related but somewhat different strategy consists in proposing an ex- 
actly solvable 2D analogue of the Luttinger model; this approach was pursued 
in [11] and [12] and again Luttinger liquid behavior up to zero temperature 
was found. 

A possible explanation of such conflicting results was suggested in [13], 
postulating the existence for the RG flow, in addition to the trivial fixed 
point associated to non interacting fermions, of a non trivial fixed point 
associated to Luttinger behavior, which could be made instable by the pres- 
ence of marginally relevant operators. In this paper we provide a quantitative 
verification of such hypothesis showing explicitely the existence of a line non 
trivial Luttinger fixed points for the RG fiow of a system of 2D interacting 
fermions with a fiat Fermi surface. It is would be not possible to derive such 
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result directly from the perturbative expansions, as it is related to cancella- 
tions between graphs to all orders of the expansion which are too complex 
to be seen explicitly; it is indeed well known that even in ID Ward Iden- 
tities (WI) are necessary to prove the existence of a Luttinger liquid fixed 
point [14]. Our analysis is based on the implementation, in an exact RG 
approach, of WI with corrections due the the cut-offs introduced in the mul- 
tiscale analysis, extending a technique already used to establish Luttinger 
liquid behavior in a large class of ID fermionic systems [15], [16] or 2D spin 
systems [17]. Such methods are the only ones which can be applied to non 
exactly solvable models, like the model analyzed in this paper. 



2 The model 

We consider a model with a square Fermi surface similar to the one considered 
in [6], [8] or [9]; the Schwinger functions are given by functional derivatives 
of the generating functional 



with i})^ are Grassmann variables, k = {k_, k^, k^), k± = j^n±, kg = 

^(?^o + 1), n±,nQ = 0,±1,±2, ... and P^di/j) is the fermionic integration 
with propagator 



_ H{k_^)Co\^a^'iki + vU\K\-pFm _ „ . . 

H{k_^) = x{a-^k-a), X{t) = 1 if t < 1 and otherwise, C^^{t) is a smooth 
compact support function = 1 for t < 1 and = for t > 7, 7 > 1. We 
assume, for definiteness, a < ao < |g- so that the support of is over 4 
disconnected regions; the Fermi surface is defined as the set of in which 
for fco = is singular, in the limit {3 00. 

By using well known properties of Grassmann integrals, see [18], ([6]) allows 
to write the Grassmann field as a sum of independent fields 

= E E (3) 

cr=it a;=it 

with "iptcryi independent Grassmann variables with propagator gu),u,v- As 
in [6], [8] or [9] we can consider only interactions between parallel patches 

'^' = EE TfTT^ E f^^(ki-k2)V^i,,,kiC2,,.,k>i,.,k3C4,a,k4^(ki-k2+k3-k4) 

a \^ h>) ki,...k4 

(4) 
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with w(x) a short range potential. The 2-point Schwinger function is given 
by 

'g2(x,y) = 1^=0 (5) 



Figure 1: The Fermi surface corresponding to the singularities of g\^] the four 
sides are labelled by (cr, u;) = (±, ±). 



3 Renormalization Group analyis 

As the interaction does not couple different a we can from now on fix cr = + 
for definiteness and forget the index a. We analyze the functional integral 
([I]) by performing a multiscale analysis, using the methods of constructive 
Quantum Field Theory (for a general introduction to such methods, see [18]). 
The propagator ([2D can be written as sum of "single slice" propagators in 
the following way 

^^(x-y)= X: e^"^-(^'+-^+)^('')(x-y) (6) 

h=—oo 

where 

and //i(k) has support in a region 0(7'^) around each flat side of the Fermi 
surface, at a distance 0(7^) from it, that is 007^^"*^ < ^0 + ^"f^X ^ ^07^''""^; 
note that in each term in ([6]) the change of variables + upp has 

been performed. The single scale propagator verify the following bound, for 
any integer M 

l9i"'(x)| < CmI^I . (8) 



X- 1 + b'^da^+l + |xo| 
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The integration is done iteratively integrating out the fields with momenta 
closer and closer to the Fermi surface, renormalizing at each step the wave 
function. After the integration of the fields ip^°\ ...,ip^^'^^^ we obtain 

I Pz,(ci^(^'^)e-^^"(^^^''') (9) 

where Pzy^idip^-^^) is the fermionic integration with propagator Z'^^{k)g^^^\ 

with g^j^^ = J2k=-oo Zh is defined iteratively starting from Zo — 1; 

moreover, if pp — {0,Pf, 0) 

OO -| 

n=l a; ^ ) ki,...,k2n i 



■ 2n 



W^ff(ki,...,k2n-l) (10) 

.1=1 

By using that / dklg^j^J^^^^] < C^^ and l^ff^^^^kl ^ C^~^ we see that the kernels 
W^^ are 0(7-'=("-2)); this means that the terms quadratic in the fields have 
positive scaling dimension and the quartic terms have vanishing scaling di- 
mension, and all the other terms have negative dimension; we have then to 
properly renormalize the terms with non-negative dimension. 

Calling k = (/c_,0,0) we define an C operator acting linearly on the 
kernels of the effective potential: 

= if n > 2 

2) If n = 1 

CW^ik) = W^{k) + kod,Mi^) + k+d+W^{k) (11) 

3) If n = 2 

CW^ik^M, ka) = Sj2^e,.„oW^{h, h, ka) • (12) 

Calling doW2{k) = —izh{k-), d+W2iy) — uiZh{kJ) (symmetry considera- 
tions are used) and 

lhik-,i,k_^2,k-,3) = (ki,k2,k3) we obtain 



= ^ U'H{k-)u;k^ - ikoZ,{k_mi^'^^^^'^ + (13) 
1 



E 7^ E ^n{k-,i, k_,2, fc-,a)v^:S^v^k1S^vJk1,53^ v^ki.5'^^(E 

w,£ I ki,...,k4 i 



where is constrained to the condition J2i ^i^iPF = and we have used 
that, by symmetry, W2{k) = 0. 
We write © as 

J p^^ (^^^{<fc))g-£v('•)(v%:^(^''')-7^v('^)(v^^(<'^)) ^^^^ 

with TZ = 1 — C The non trivial action of TZ on the kernel with n = 2 can 
be written as 

Wi^(ki, k2, ka) = [W^{}ii.. k2, kg) - W^{}iiM, ks)] (15) 

+ [W^{\^^MM) - W^4'(kl,k2,k3)] + [Wi^(ki,k2,k3) - W^i^(ki,k2,k3)] 

The first addend can be written as 

^1 / (it4(,iiy4(A;_,i,/c+,i,tA;o,i;k2,k3) + 
Jo 

k+,1 /'rft9fc+,,W^i'(^-,i,^^+,i,0;k2,k3) (16) 



The factors fco,i and k+^i are 0(7^'), for the compact support properties 
of the propagator associated to 'ip^^^\ with h' < h, while the derivatives 
are dimensionally 0(7"'^"^); hence the effect of TZ is to produce a factor 
7^ -'i-i ^ I making its scaling dimension negative. Similar considerations 
can be done for the action of TZ on the n = 1 terms. The effect of the C 
operation is to replace in 14^2'' (k) the momentum k with its projection on the 
closest fiat side of the Fermi surface. Hence the fact that the propagator 
is singular over an extended region (the Fermi surface) and not simply in 
a point has the effect that the renormalization point cannot be fixed but it 
must be left moving on the Fermi surface. 

In order to integrate the field ip^'^^ we can write (fT4l) as 

J Pz, ^((i?/>(^^))e~'^'^''^^'^'^"')"^^"''(^'^*^"') (17) 

where Pz^_^{dil)^-^^) is the fermionic integration with propagator 

1 H{k_)C^\\^) 
Z/i_i(k) —ik^ + uJVpk^ 

with C,-i(k)=Et-oo/fe and 

Z,_i(k) = Z^(A;_)[1 + Hik_)C^^\k)zhik_)] (19) 

Moreover CV^ is the second term in (|T3l) . 
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We rescale the fields by rewriting the r.h.s. of ffT^ as 



where 

_ * 1 _ _ _ _ 

^ I ki,..,k4 j 

(21) 

and the effective couplings 



^"^-"^ Uk^,i,k_,„k_,,) (22) 



^=1^ Zh-iik^^i) 

After the integrations of the fields 'ip^'^\ ^p^~^\ -0^'*^ we get an effective the- 
ory describing fermions with wave function renormalization and effective 
interaction fl^Tl) . Note that and gh are non trivial functions of the mo- 
mentum parallel to the Fermi surface. 
We write 



Pz,,,(#^-'-')) j P^^_^(rfV^W)e-^^''"(V^^''"'')-^^'''(V^^<-'^') (23) 
and the propagator of Pzi^_-^{d'ilj) is 



Zh-i{k-) —iko + LUVpkj^ 
and 

with H{k_)fh(k) having the same support that H{k^)fhik). We integrate 
then the field ijj^'^'^ and we get 



Pz,_,(^^V'(^'-'))e-^"'"''(V^^*''"') (25) 

and the procedure can be iterated. 

The above procedure allows us to write as a series in the effective 
couplings gk, k > h, which is convergent, see [18], provided that L_ is finite 
and Eh = s^Pk>h\\9k\\ small enough; moreover ||W2n|| = 0(7"''^""^^). A 
similar analysis can be repeated for the 2-point function. 
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However even if the couplings gk starts with small values, they can possi- 
bly increase iterating the RG and at the end reach the boundary of the (es- 
timated) convergence domain; if this happen, all the above procedure looses 
its consistency. A finite temperature acts as an infrared cut-off saying that 
the RG has to be iterated up to a maximum scale hp = 0(log/3) and, up to 
exponentially small temperatures i.e. P < 0(e'^''^' ), then surely the effec- 
tive couplings are in the convergence domain; however, in order to get lower 
temperatures, more information on the effective couplings are necessary. 



4 RG flow and Luttinger liquid flxed point 

The RG analysis seen in the previous section implies that the effective cou- 
pling Qh verify a flow equation of the form 

gh^i = gh + Pi''\gh;-;go) (26) 

where the r.h.s. of the above equation is called Beta function, which is 
expressed by a convergent expansion in the couplings if Eh is small enough. 
The first non trivial contribution to /3g^\ called P^^^^'^\ is quadratic in the 
couplings and it is given by 

/3f ('^^ = + (27) 

where 

PI^^ = J dpH{h,. ~ p^)H{h,- + p-)gh{h,.,h,. -v-.h,-) (28) 

(1 11., Jhip)Chip) 

9h[ki- -P-,k2-,h.- +P-)^-, — 

Po + VfP+ 

Ph^ = - J dpH{k2- -p^)H{ks_ - P-)ghiki-,k2- -P-,k3,- - P-) 

9h{k2,~ - A;2,-, A;3,-)^- — (29) 

p^o + vppi 

The above expression essentially coincides with the one found in [6] or [8] ; it 
is indeed well known that the lowest order contributions to the Beta function 
are essentially independent by RG procedure one follows. 

The flow equation (!26|) encodes most of the physical properties of the 
model, but its analysis is extremely complex. Some insights can be obtained 
by truncating the beta function at second order, and by the numerical anal- 
ysis of the resulting flow by discretization of the Fermi surface; it is found, 
see [6] or [8],that gh{k^^i, k^^2, k-^s) has a flow which, for certain values of 
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A;_ 1, A;__25 ^-,3 increases and reach the estimated domain of convergence of 
the series for W^''\ While this increasing can be interpreted as a sign of in- 
stability, mathematically speaking this means that the truncation procedure 
becomes inconsistent. 

A basic question is about the fixed points of the flow equation ( l26l) : in 
particular if there is, in addition to the trivial fixed point Qh = 0, a. non trivial 
fixed point corresponding to Luttinger liquid behavior. Note first that the 
set ^ 

gh{k-A, A;_,2, k_^s) = —S{ki_ - k2-)\h (30) 

with \h constant in k, is invariant under the RG flow, in the sense that if 
£V'''^^ has the form,for k > h 

^/»<^p-,o^<::^,k^j,k+pV'<i'',k'^j',k'-p (31) 



k,k',p 



the same is true for £V(^~^). This can be checked by the graph expansion. In 
the graphs contributing to wl'^^ , the external lines of the graphs contributing 
to 1^4'*^ either comes out from a single point, or are connected by a chain of 
propagators with the same a;, Moreover in each Feynman graph the only 
dependence from the momenta of the external lines is through the function 
H{kJ) which are 1 in the support of the external fields / dkH(k-)'ip^ = 
Jdkij^. For the same reasons also Zh is independent from A;_. 

The crucial point is that, in the invariant set (!30|) . some dramatic can- 
cellation are present implying the following asymptotic vanishing of the beta 
function (which will be proved in the subsequent sections) 

/3f ) = Oi^'^el) (32) 

saying that there is a cancellations between the graphs with four external 
lines and the graphs with two lines contributing to the square of Zh , see 
( I22I) : such graphs are 0(1) but there are cancellations making the size of the 
sum of them 0(7'^). At the second order (1321) can be verified from fl29p and 
(IHUj): at third order, it is compatible with (A16),(A15) and (4.8) of [8]. 

The validity of ( |32l) immediately implies the existence of a line of non- 
trivial fixed points for (l26i) of the form 



g^oc{k-,l, fc-,3) = -^(5(fci _ - /C2 -)A_oo (33) 

with A_oo k-independent and continuous function of U, A_oo = Aq + 0{U'^), 
Xq = cU for a suitable constant c. 
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Note also that to such fixed point is associated Luttinger hquid behavior, 
as, from (fTUl) . ^ 7^'''', with 77 = cAq + 0{U^) and p = (0,p+,po) 

a = hm l^YYT dp+dpo — ; — — 34 

Hik'_)aik' + p) d H{k_)Cj;\k-p) 

-i{kQ + po) + uj'vF{k+ + p+) dk+ ~i{ko - po) + uJVF{k+ - p+) ^=0-^+-° 
Indeed the 2-point Schwinger function can be written as 

<^=± k 
with AW(k) = 0{eh), so that 

,S2(x,y)= ^ e^-pH-+-?;+) (35) 

^ ^ i/(A:,)Co-^(k) 1 + A(k) 



V(3 Y -^^0 + 0;^;^/?+ \kl + vpil'' 

with |v4(k)| < C|[/|. This means that to the fixed point is associated Lut- 
tinger hquid behavior, as the wave function renormahzation vanishes at the 
Fermi surface as a power hke with a non-universal critical index;the Luttinger 
liquid behavior is found only if L_ is finite, as if L_ — >^ cx) the critical index 
is vanishing. Note also that the cancellation in ( l32l) reduce to the one in ID 
if /e_ = k' in (EB. 



5 The auxiliary model 

There is essentially no hope of proving a property like fl32p directly from the 
graph expansion, as the algebra of the graphs is too cumbersome (except 
than at one loop in which it is easy to check). We will follow instead the 
same strategy for proving the asymptotic vanishing of the Beta function in 
ID followed in [15], [16], considering an auxiliary mode/ with the same beta 
function, up to irrelevant terms, but verifying extra symmetries, from which 
a set of Ward Identities can be derived. In the present case, such identities 
are related to the invariance under local phase transformations depending 
on the Fermi surface side, which in the model ^ is broken by the cut-off 
function Co'^(k) and by the lattice. 

We consider an auxiliary model with generating function is 

(36) 
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with Q-jv(k) = Eth fk{kl + vlkl), h < p^,. = tA^^^x,. 
- U r 

= Y~ (RV\'i ^ lrr~ ^ ^(^^0,P+)^t!^,k^c^,k+(0,p+,po)'^'^^k'''/^L..',k'-{0,p+,po) (^'^) 

^- I k,k' '^-^+ po,p+ 

with v{xq, x+) a short range interaction. The above functional integral is very 
similar to the previous one, with the difference that there is an ultraviolet 
cut-off 7^ on the + variables , which will be removed at the end, and an 
ultraviolet cut-off 0(1) on the — variables; such features are present also in 
the models introduced in [9] or [12]. 

Again (l36il can be analyzed by a multiscale integration based on a de- 
composition similar to (E]), with the difference that the scale are from h to 

In the integration of the scales between and 0, the ultraviolet scales, 
there is no need of renormalization; apparently the terms with two or four 
external lines have positive or vanishing dimension but one can use the non 
locality of the interaction to improve their scaling dimension. We integrate 
(with £ = 0) the fields and we call M^£!m the kernels in 

the effective potential multiplying 2n fermionic fields and m J fields. Again 
the dimension is 'y-'=("+™-2)^ A; > and we have to improve the bounds using 
the non-locality of the interaction. We can write 

<(x, y) = / rfyi^ ^^"°""^';^'_""~"^'"V ^?(yO.^^-")(x - y.)<(y2; y) 

+U j dy /^"° - y^'l - '-''^ ,(^-^)(x - y,)M^g(y, y,; y,) + (38) 
f/5(x-y) /dy /^"°-^^-f^-^^-^V o^?(yO 



The first and the third addend of Fig. 2 are vanishing, by the symmetry 
5r(fco, /c+, kJ) = —g{—kQ, —k^, k_); hence, using that < Cj'^ and that 

1^2^'' is 0{U) (by induction), we obtain the following bound 

IliyiSlI < • E \\g^% < C\U\^'L-Jr'' (39) 

j=k 

Note that we have a gain 0(Ll^7~^'^)), due to the fact that we are integrating 
over a fermionic instead than over a bosonic line. 
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Figure 3: Decomposition of W2J 

Similar arguments can be repeated for ^ which can be decomposed 
as in Fig 3. The second term in the figure is bounded by 0(|t/|Ll^7~^'^). A 
similar bound is found for the third term in Fig.3; regarding the first term, 
we can rewrite it as 

J d^dzlg^'^'^Xz - ^)f^v{x, - zo, x+ - y) = 

/ rfxdz^ "^"'° - - [,(^^^)(x - z)]^<j(z, y) (40) 

+ / d:sidz—[v{zo - xq, z+ - x+) - v{zo - zq, z+ - z+)][g^^'^\y. - z)fW^^i{'z, yo) 
and using that 

J dx[^^('='^)(x - z)]2 = J dk^H{k^) J dk,dKj^^^-y^ = (41) 

the first addend is vanishing; the second addend, by using the interpolation 
formula for v{zq — Xo,^+ — x+) — v{zq — zq,z+ — z^), can be bounded by 
C\U\'y~'^, as by induction ||Wo^2^|| < C\U\. A similar analysis proves the 
bound for W^g. 
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After the integration of the fields tp^^\'ijj^^~^\ ...,tp^~^^ we get a Grass- 
man integral very similar to ([T]); the integration of the remaining fields 
:0{o)^^(-i)^ .. is done following the same procedure as in section 3, with the 
effective coupling of the form ( l30i) and ^V*^*^-* of the form ( l3Ti) . The crucial 
point is that the beta function coincides with the beta function for the model 
([T]) up to 0(7'^) terms; hence it is enough to prove the validity of in the 
auxiliary model. 

6 Ward Identities 

We derive now a set of Ward Identities relating the Schwinger functions of 
the auxiliary model ( l36l) : by performing the change of variables 

<x ^ (42) 

and making a derivative with respect to ax,^; and to the external fields we 
obtain 

J dk'[H-\k'_ +p_)ChA'^' + p){-i{k',+po)+u;vF{k'^+p+)) - 

H-\k'_)CnA^){-ik', + ^t;i.fc:,)](V^+^,+pC,k'V^.^',k-pC',k) = 
5.,.'[(^:^',k-pfc,k-p) - {^^',^3Z',^.)] (43) 

where (^:i",k'+pC,k'^^',k-pC',k) is the derivative with respect to Jp, 0+,,k_p, C',k 
of flHUl) . Computing fHSj) for p_ = we get, if p = (0,p+,po) 

{-m + ^^FP+)(pp,<^V'k,<^'^k„p,^/) + A(k, p) = 

[(^ci'',k-p^j',k-p) - (v^:^',kC',k)] (44) 

and 

A(k,p) = j c^k'C(k',p)(VJ+,,^pC,k'^:^',k-pfc,k) (45) 

with 

C(k, p) = (-ifco + uvFk+) (46) 
[Cfe,^(k + p) - C/,,7v(k)] + (-ipo + c^^^fP+) [ChA'^ + P) - 1] 

In deriving the above equation we have used that 

/ rfk'i/-i(A:'_)(ft_,,^p^-^,V^+,i,_pV^-,) = 
/ rfk'(^^+k'+pfck'^:r',k-pfc,k) (47) 
13 



for the compact support properties of the fields ?/^k and = H; note also 
the crucial role of the condition p_ = in the above derivation. 

The presence of the term A(k, p) in the Ward Identity (jHj) is related to 
the presence of the ultraviolet cut-off; as in ID, such a term is not vanishing 
even in the limit N ^ oo and it is responsible of the anomalies, see [15], [16]. 
The following correction identity holds, similar to the one in the ID case 

A(k, p) = u{-ipo - UVFP+) E (Pp,-"V^:i'',k-pC',k) + (k, p) (48) 

w"=± 

with a small correction. Indeed can be written as functional deriva- 
five, with respect to 0"*", 0", J, of 

gWA(J,0) ^ j p{^dilj)e~^^'^^^^'-^ rfz['/'i,z</'J,z+<zV'J,z]+To(J,V')-T_{j,v) (^49^ 

with 

T-m -Ej ji^)^)j^,-'M-ip« - -'»'FP+)5ir+p,„.fe. (51) 

(H9l) can be evaluated by a multiscale integration similar to the previous one, 
the only difference being that / Jp is replaced by Tq — T_. The terms with 
vanishing scaling dimension of the form Ji/j^i/j' can be obtained from the 
contraction of Tq and T_; in the first case we can perform a decomposition 
similar to the one in Fig.3, see Fig. 4. Regarding the second and^^hird Jerm, 




Figure 4: Terms obtained from the contraction of Tq; the black dot represents 
C(k,p) 

we can proceed exactly as in the previous section, the main difference being 
that at least one of the two fields have scale so that they obey to the bound 
0('y-2fc^-(i/2)(Ar-fc)'j_ 'jj-^jg foUows from the fact that when C is multiplied by 



14 



two propagators we get 

/.(k+P) ,/.(k) 



which is non vanishing only if one among i or j are equal to h or A^. 

The main difference with the analysis in the previous section is that in 
the first term of Fig 4; the "bubble" in Fig. 3 was vanishing, while here it is 
not. We choose in (HUI) equal to the value of this bubble, in order to cancel 
it. The value of the bubble is given by 

- = Uvip„p^) f ^\ _ ^i^^Hkki^^Hk + p) (53) 

J {27V)'^ —ipo — UJVFP+ 

and, in the limit N oo, ly = Uv{j)Q,p+)^. Hence for fco, A;+ = 0(7'') 

\R'^\\.,t>)\<Ceh-'' (54) 



7 Schwinger-Dyson equation 

An immediate consequence of the analysis in the previous section is that, for 
momenta computed at the infrared scale |k| = |k'| = |k + p| = |k' — p| = 



1 ih) (h) ih) (h) ^( XX 

— ^c.,k^c.',k'^^,k'+5^^',k'-5<>^ + '-'y^h)) 



Ah) 

(C,kfc) = ^(1 + 0(6.)) (55) 

This says that relations between the effective couplings at a certain scale 
h < Q can be obtained from relations between the Schwinger functions of 
the auxiliary model ( l36l) computed at the infrared cut-off scale. The starting 
point for deriving such relations is the Schwinger-Dyson equation for the 
4-point function, given by, if the external momenta are computed at the 
infrared scale and p = (0,p+,po) and p = (0,p+,po) 

(^j,k^:^,k-5fc,k'^^,,k'+5) = 

^{— t;(po,P+)£/^^k'+5(C',k'^<:^',k')(P5,'^''fckV'^,k-5) + (^6) 

U)" ~ 

/ ^^^(^'o>l'+)(Pp,^''fckV^:^,k-5fc,k'^^^k'+5-p)} 
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Figure 5: Graphical representation of fl56|l : the dotted hne represent the free 
propagator 

By the WI dH,® 

(P) - (Ck^^,k)] + (k, p) (57) 

with A^^^{p) = 1 + 0{£h), A^-uj{p) = 0{eh) and even in p; moreover 
-ffa;,w'(k, p) is a hnear combination of the i?^'^ functions in ( l48l) with bounded 
coefficients. The WI for the 4 point function is given by 

UJ' 

and similar ones, so that the second addend of the l.h.s. is given by 

J dPodp+XeiP) (^k-p,..^k-5,..^k',..' V^k.+5_p,..) 

+MP){kJt^J^',^'^t+^,J + ^'''(k,k',p)] (59) 

where Xe{p) is a compact support function vanishing for p = 0, and such 
that it becomes the identity in the limit e — 0; moreover the functions Ai{p) 
are bounded and even in p, and i?^'^ is a linear combination of the i?^'^ 
functions in (!58|) with bounded coefficients. 
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We have now to bound all the sums in the r.h.s. of fl3U|) . Note first that, 
by parity 



/ dPodp^Xe{P)AM-l^f^;^ = (60) 

Moreover the first term in the r.h.s. of (|59ll verifies 

(61) 

and a similar bound is true for the second and third term. Finally, as in [16] 

I / dpodp^ ^iPo,P+) H'A^k,k',p)\ < Ce,^ (62) 
J -ipo + vfuj"p+ 

By inserting ([55]), (IS?!), (El), ([60]), (E^ in §B we get = Aq + 0{U'^), which 
means that the effective interaction remain close to initial value for any RG 
iteration; a contradiction argument shows that this can be true only if the 
beta function is asymptotically vanishing: as this beta function is the same 
of the 2D model ([T]) with effective couplings (30), up to 0{'y^) terms, then 
(|321) follows. 



8 Conclusions 

We have shown that the RG flow for a system of spinless fermions with flat 
Fermi surface has, in addition to the trivial fixed point, a line of Luttinger 
liquid fixed points, corresponding to vanishing wave function renormalization 
and anomalous exponents in the 2-point function; such fixed point is in the 
invariant set (30). This makes quantitative the analysis in [10], in which 
the existence of a Luttinger fixed point in 2D was postulated on the basis 
of bosonization. With respect to previous perturbative RG analysis,the key 
novelty is the implementation of WI at each RG iteration, in analogy to what 
is done in ID. 

Of course the other effective interactions should cause flows away from 
this fixed point. Indeed the situation is somewhat similar to the ID (spin- 
ning) Hubbard model, in which there is a Luttinger liquid fixed point in the 
invariant set obtained setting all but the backscattering and umklapp scatter- 
ing terms equal to zero (that is the set gi^h = g3,h = in the g-ology notation, 
see [19]), and a flow to strong coupling regime driven by the backscattering 
interaction. We can in any case expect, as in ID, that even if the Luttinger 
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fixed point in 2D is not stable its presence has an important role in the 
physical properties of the system. 

Fermi surfaces with fiat or almost flat pieces and no van Hove singularities 
are found in the Hubbard model with next to nearest neighbor interactions 
or in the Hubbard model close to half filling, and it is likely that our results 
can be extended, at least partially, to such models. Note however that in 
such models the sides of the Fermi surface are not perfectly flat, so that 
one expects a renormahzation of the shape of the Fermi surface, as in [3], 
which is absent in the case of flat sides by symmetry. Another simplifying 
property of the model considered here is that the modulation of the Fermi 
velocity is taken constant along the Fermi surface, contrary to what happens 
in more realistic models; a momentum dependent Fermi velocity produces 
extra terms in the WI, as it is evident from (43) {vp should be replaced by 
vpik'J) in the first line and VF{k'_ -\- p-) in the second line), and their effect 
deserves further analysis. 
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institute in Wien where this paper was partly written. 
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